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1 Introduction 



A routine numerical solution of an asymmetric Schrodinger bound-state problem on 
the line x G (—00,00) requires a careful verification [|l|. One needs non-numerical 
asymmetric models. For this purpose we may use the shifted harmonic oscillator, 
Morse's well and the two scarf-shaped hyperbolic forces. All of these models (cf. 
Table |l|) are listed in review as possessing the complete solution in closed form. 

There exist incompletely solvable polynomials V{x) = ax + bx'^ + ... + zx^ 0] 
and multi-exponentials V{x) = ae~^ + 6e~^^ + . . . + ze~^^ [Q. They extend the 
possible tests and further non-numerical applications beyond N = 2. In a puzzling 
contrast, a natural generalization 



Mr M 
m=Z n=i 

of the remaining two items in Table |l| is not amenable to the similar elementary 
treatment |^ . This distracts attention from the hyperbolic oscillators ([1|) in spite of 
their obvious phenomenological as well as purely mathematical appeal. 

In the present paper we shall return to several formal as well as descriptive par- 
allels between the separate items in Table |I[ On their basis we shall propose and 
describe a new semi-analytic approach to the "neglected" family (|^). 

In Section 2 we recall the harmonic and Morse oscillators and their N > 2 gener- 
alizations as our overall methodical guide. In the language of the well known Lanczos 
method we underline the key role of simplicity of the repeated action of the Hamil- 
tonian upon a suitable trial state |0). An appropriate choice of this initial ket vector 
is able to inspire some of the existing non-numerical power series solutions. In this 
setting the Lanczos approach is shown to find its natural re-incarnations in the well 
known method of Hill determinants [0] as well as in the symmetric Jost-solution 
method of ref. p[. 

In Section 3 we show that and how the latter two examples pave the way towards 
eq. ([l|) with any M > 2. In a full parallel to the polynomial case we construct the 
asymptoticaly correct bound state solutions which all retain a recurrently defined 
power-series structure. Via an appropriate D— dimensional partitioning of the basis 
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we preserve their connection to the two remaining exactly solvable hyperbolic M = 2 
examples of Table |T]. 

Section 4 illustrates the technical details at the first nontrivial D = 2. We 
contemplate there a spatially anti-symmetric M = 2 exercise (|I]) using /2 = S'l = 0. 
We detail the proof of the point- wise convergence of our "partitioned hypergeometric" 
wave functions. We show how the symmetry considerations significantly simplify the 
construction and matching of our wave functions near the origin. 

Section 5 adds a short summary. 



2 The method 

2.1 Wave functions in the Lanczos basis 

The Lanczos numerical eigenvalue method works with a set { |n) } of the basis 
ket vectors which are generated via a repeated action of the Hamiltonian H upon 
an initial vector |0). In a slight generalization of this procedure one has to assume 
that the action of the full Schrodinger operator H — z upon each ket \n) may be 
represented as a linear superposition over the same set of the kets 0, 



[H - z)\n) = |0) ■ Qo,nW + |1) ■ Ql,nW + 



With a matrix of functions Qm,n{z) (cf. p. 257) we may abbreviate 



(2) 



iH-z)\0),iH-z)\l),... 



|0),|1),|2),... 



(H-z) 



X 



(H-z) 



X 



X 



and solve any linear homogeneous equation {H — E)\y) = hj the ansatz 

CO 

\y) = Y: \n) K ^\X}h. (3) 



n.=0 



Provided that the separate lanczosean kets are linearly independent the resulting 
identity \X}Q{z)h = may be interpreted as a system of conditions 



Q{z)h = 



(4) 
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The practical applicability of this recipe relies upon several tacit assumptions. Most 
often one chooses the set { |^) } as a common harmonic oscillator basis ||TT|. It is 



orthonormal ({X|X} = /) and complete {\X} {X\ = Id) and we may truncate the 
linear set (^) to the mere routine matrix diagonalization 

M 

T.[Qi^)-ElUnhn = 0, m = 0,l,...,M, (5) 

n=0 

This is a textbook variational recipe and its secular equation 

detQ{E)=0 (6) 



determines the spectrum numerically |^ 



A non-variational and less numerical modification of the construction may be 
based on a more sophisticated choice of the Lanczos basis. Various linear algebraic 
algorithms of such a type are used to solve various Schrodinger equations in appli- 



cations [13|. Let us recall two examples as our methodical guide. 



2.2 Anharmonic example 

Both the above-mentioned multi-exponential and polynomial oscillators prove mutu- 



ally equivalent after a change of variables [0]. Their "canonical" [15[ representation 



^(^) = 7^ + 9ir' + 92r^ + ■■■+ g2N-ir^''-\ r E (0, oo) (7) 
is easily tractable by the variational algorithms. In the less numerical power-series 



approaches [T^ the harmonic kets are being replaced by their mere power-law com- 
ponents {r\n) = (r|0['^'^l) ■ r". This leads to an asymmetric matrix Q{z). Its linear 
algebraic eq. (^) proves often solvable as a very simple recurrent specification of the 
coefficients hn in eq. (|^) (cf. ref. [|l^ for more details). 

An even more ambitious reduction of Q may be achieved after an anharmonic 
choice of the initial |0). According to Magyari this assigns a few elementary 
bound-state solutions to many multi-exponential and polynomial potentials at cer- 
tain exceptional couplings. At arbitrary couplings and energies the same option |0) 
may provide an extremely compact infinite-dimensional algebraic secular equation 
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(^). For illustration let us consider the famous sextic oscillator example of ref. |18|. 
With = 2 in eq. (|^, denoting = 16 and g2 = 16 a P and using the WKB- 
inspired postulate 

{r\n) = ^n+l+l^-c^r^-pr^ ^ ^ ^ q 



(8) 



we get the tridiagonal quasi-Hamiltonian 

/ 

Q{E) = 

V 





71 












72 










^2 


73 



(9) 



Its equation (|[) may safely be interpreted as an infinite-dimensional limit of the 
truncated diagonalization provided only that (72 > ||19[. The three nonzero 
diagonals in eq. (§) have to be compared with the seven-diagonal structure of the 
Hamiltonian in the usual orthogonalized harmonic oscillator basis. 

For (72 ^ and at a special discrete set of the couplings gi the infinite-dimensional 
tridiagonal secular Hill determinant factorizes and the recipe reproduces a part of 
the spectrum correctly ||18|. In all the other cases the WKB-compatible Lanczos 
basis ceases to be adequate. The Hill-determinant recipe loses its relation to the 
correct asymptotic boundary conditions and the basis (0) must be regularized for 



certain hidden-symmetry reasons [^. More diagonals necessarily appear in eq. (^. 
Otherwise, one gets wrong results from the truncated eq. even in its infinite- 



dimensional limit W% - 

Virtually no similar constructions of our short-range hyperbolic oscillators seem 
to appear in the current literature. In the present paper we intend to explain the 
difference and develop a new semi-analytic approach to eq. (|l|). Our construction 
will fairly closely parallel the formalism of the Hill-determinant method. In our 
second preparatory step the appropriately modified choice of the Lanczos basis will 
be illustrated via the symmetrized Rosen- Morse or scarf model of Table |l[ 
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2.3 Poschl- Teller example 



Formula (|l|) with M = 2, attraction / = — A(A — 1) and vanishing g = defines the 
bell-shaped and spatially symmetric Poschl- Teller well V^-^'^\x) = //cosh^x p2 |. 
The functional form of the optimal lanczosean kets is more or less uniquely deduced, 
very much in the spirit of the "most ambitious" WKB-like choice in eq. (H) above, 
from the available exact solutions, 

sinh^^-? 



X 



G L2( — oo, oo) . 



(10) 



cosh-^+^^+Pa; 

All these basis states possess the even or odd parity at g = 1 or g = 0, respectively. 
Within this subsection let us fix p = 1 — g. Then, the action of the full Hamiltonian 
ff{PT) _ _g2 _j_ v^^'^\x) on our symmetrized/anti-symmetrized states ( p!OD becomes 
particularly transparent. For energies E = —k"^, it is characterized by the mere 
two-diagonal matrix 

/ 



Q{E) 



ao 

/3i a2 









V 



with the vanishing uppermost element ao = 0. The bound-state solutions @ of our 
Schrodinger differential equation read 

As long as they are defined by the really elementary two-term recurrences (j^) 



riQ 



|0)-|1).^ 



(12) 



/ 









-4(k + 1) 
f + {k + p + 2){k + p + 3), -8(k + 2) 



V 



hi 
h2 



\ ■ J 

(13) 



our solution \y) coincides with the Gauss hypergeometric series. 



hn tanh^ x 



1 

cosh'' X 



K+p+X K+p+l—\ 



1 + k; 



cosh X 



(14) 
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It is defined on a half-axis, say, x > 0. Fortunately, due to the manifest symmetry or 
anti-symmetry of the physical solutions the necessary analytic continuation across 
the origin proves equivalent to the termination of this infinite series. The well known 
Jacobi polynomial solutions are obtained at each physical energy p[. 

3 Partitioned expansions 

We may conclude that the description of bound states by the infinite series (U) proves 
easy and efficient not only in the Hill-determinant setting of section 2.2 but also in 
an alternative Jost-solution spirit of section 2.3. We intend to extend the parallelism 
far beyond the trivial example of section 2.3. 

The action of the kinetic energy T = —d1 on the basis ( p!OD conserves both the 
independent parity-like parameters p and q. The same conservation law is obeyed 
by the single-term symmetric potentials V^'^^\x) = f /cosh'^^ x with the even ex- 
ponents M = 2K. The rule is broken by the general Hamiltonians containing 
superpositions of the symmetric and ant i- symmetric components Vj'^^\x) and 
V^^\x) = sinh x/cosh^ X, respectively. Nevertheless, the full basis ( p!0D numbered 
by a composite index /i = fi{n,p,q) = An + 2p + q > 1 (as in,p,q,K{.x) = 
/i = 1, 2, . . .) proves reducible for all the single-term potentials V^^\x) = ±.V^^\—x) 
of a definite parity. 

3.1 Symmetric potentials V{x) = V{—x) 

We may choose the initial Lanczos ket |0) either as the spatially symmetric (and 
asymptotically correct) element (x|S^(o,o,i)) = cosh'^'x with p = and g = 1 or as 
its anti-symmetric analogue (x|S^(o,i,o)) = sinhx ■ cosh~'^^"'^x with p = 1 and g = 0. 
In both these cases, all the Hamiltonian operators T -|- V"/^'^-* (x) become compatible 
with recurrences (|^) in the two alternative bases 

|0)) |2), • • • = |2^(o,o,i)), |2^(i,o,i)), |2^(2,o,i)) • • • = IS5), IS9), . . . , 

10)5 11)5 . . . = 1^^(0,1,0))) |'='/i(2,l,0)) • • • = |'='2), |'='6)) |'='lo); • • • 
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with p— 1 — q — OotI, respectively. After we abbreviate Oj = —j {2k + j) and 
bj = {k + j){K + j + 1), this enables us to reproduce the two-diagonal Poschl-Teller 
realization oi Q = Q^^^ at = 1, 

^ ... 

f + bo, 02 ... 

/ + &2, 04 







V 



/ 

f + bu 



I 









02 
/ + &3, 04 



V 



In the "first unsolvable" case with K — 2 the couphng / moves one step down, 



/ 



V 












... 









... 


/ 


&2, 


04 


... 





/ 


64, 


06 



/ 



Q 



(1) 



V 












... 




02 





... 


/ 




04 


... 





/ 


^5, 


06 



Partitioning indicated by the auxiliary lines tries to preserve the same two-diagonal 
pattern as above. At = 3 we have, similarly, 

/ 



Q 



(0) 















\ 




/ 














\ 




02 














bi, 


02 













h2, 


04 














bs, 


04 







/ 





64, 


06 








/ 





bb, 


06 








'■■ ) 




\ 







V 

and so on. The dimension of partitions grows linearly with M = 2K as D = K. 

The second series of the symmetric potentials V}^^^~^^\x) — f/ cosh^^"'"^a; with 
the odd powers M — 2K + 1 must be investigated separately. It acts on our parity- 
preserving basis in such a way that the conservation of the quantum number p is 
broken. The following g— preserving bases must be used. 



|2), . 


• • = ■^2), ■=^4), l^e), l^s), ■ ■ ■ , 


g = 0. 


|2),. 


■ ■ = Si), S3), S5), S7), . . . , 
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At = the new lower triangular matrices Q = Q'^' contain just the three nonzero 
neighboring diagonals. For a preservation of the two-diagonal denotation it is suffi- 
cient to switch to the D = 2 partitioning. Similarly, a three-dimensional partitioning 
is needed at i^' = 1. With the further increase of K the dimension D = 2K+1 grows 
more quickly. 

3.2 Anti-symmetric potentials V{x) = —V{—x) 

The class of the anti-symmetric forces V^'^^\x) with L > 1 inter-relates the basis 
states with different parities q. The value of the index p is conserved, 

Vl'^^\^)\'^Kn,P,q)) = (1 - fl' ■ \'^fi{n+L-l,p,l-q)) + " | S^(n+L,p,l-g) ) • 

The Hamiltonian T + V^^^ acts transitively on the following two reduced Lanczos 
bases, 

|0), |2), . . . = IS4), IS5); iHg), IS9); . . . , p = (15) 
|0), . . . = IS2), IS3); iSe), IH7); |Hio), |Sn); . . . , p = 1. (16) 

Marginally, we may note that at L = the structure of the matrix Q ceases to be 
triangular. This seems closely related to the asymptotic asymmetry of the g\ ^ ^ 
potentials V[^'^l(— 00) = —gi ^ V"[*^l(oo) = +gi and to their anomalous non-Jost 
solvability via a change of variables at M = 2 (cf., e.g., [^). In this subsection 
we shall assume that gi = 0, therefore. This constraint is further supported by the 
observation that at M = 1 the monotonic Vj-"*^^ {x) = gi ■ tanh x itself cannot generate 
any bound states at all. Thus, our study of the anti-symmetric models has to start 
at the exactly solvable Vlf'\x) = (7 sinh x/cosh^ x (cf. Table p. 

This anti-symmetric scarf (AS) potential V^'^\x) = V^'^^^^x) is extremely suitable 
for methodical purposes. Its significance is connected to the fact that our basis ([10|) 
is not tailored precisely to its exact solvability. A D = 2 partitioning is needed. 
In the reduced bases (^) and (|r^) its recommended boundaries are marked by the 
semi-colons. For all the L = 2, 3, . . . descendants V^'^^\x) of the AS example the size 
D of partitions will grow due to the downward shift of the constant g again. 
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The action of the last class V^'^^^^''{x) = g ■ sinhx ■ cosh~^^~^ x of the simplified 
single-term potentials on the kets (|T^) looks irreducible. The impression is wrong. 
After we introduce a new quantum number I = 2p + q (modulo 4), the basis elements 
with 1 = and 1 = 3 never mix with their 1 = 1 and 1 = 2 counterparts. For both 
the initial choices of |0) = and |0) = IH2) we arrive at the same output, 

10)5 11)5 |2), ■ ■ ■ = or 2), 1^5), l^e), 1^9), \^lo), l^is), • • • • 

The difference between the two matrices Q will only lie in their elements. 
3.3 Asymmetric Lanczos kets 

Asymmetric oscillators (0) admit a non-conservation of parity by each Lanczos ele- 
ment \n) separately. The functions 

{X\n) = e„,p,,,a,«(x) = ^l^llZ^ arctan(.nh.) ^ ^) (^7) 

generalize their a = predecessors (p!OD and represent a very good new candidate 
since, due to the presence of a new parameter a, the number of the new terms in 
eq. d^) may be lowered, for any potential (|l|), more efficiently. First of all, this 
implies that we may admit the nonzero gi again. Via a suitable choice of the value 
of a we shall be able to reproduce all the "missing" (viz., Rosen Morse and scarf) 
terminating solutions of ref. or Table [l|. 

At a 7^ also the action of an arbitrary hyperbolic Hamiltonian remains trans- 
parent and elementary in the purely kinetic limit, 

^nnoo«-(^) / n9 — cr(cr + 1) — (2(7 + l)a sinh X , .o — 2asinhx 
= (a + g-lf + ^ ^ + (g-l)^— -2 . 

C,n,p,q,a,K{^) COSh X Smh X 

Here, a = a{n,p) = K, + 2n + p and the prime denotes the differentiation with respect 
to X. The action of the purely kinetic Hamiltonian T = —d'^ on our innovated kets 
= C,n,p,q,a,K{x) may employ the multi-indices fi{n,p, q) = in + 2p + q again, 

T\E^(n,p,0)) = — (O" — 1)^ |S^(n,p,0)) + (2CT — 1) fl | -|- 

+ {a'^ + a - a^) |E:^(„+i,p,o)) - (2cr + 1) a |E:^(„+i,p,i)), 
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^|S^(n,p,l)) = — O"^ + (2(T + 1) a |S^(„+i^p^o)) + ((T^ + Cr — a^) \'^fj,(n+l,p,l)) ■ 

The kinetic matrix elements of Q depend on a and a and all of them increase with 
n. Due to the presence of the new parameter a the kinetic operator T inter-twins 
the states (0) with different parities g = 0, 1. The states with different p = 0, 1 
stay decoupled. 



3.4 Partitioned hypergeometric-like series 



Our present proposal may be summarized as an application of expansions @ to 
potentials (|1|) inspired by the analogies between the Poschl- Teller and harmonic 
oscillators. The feasibility of our construction stems from the fact that the action of 
the present class of Hamiltonians on the suitable Lanczos kets may be characterized 
by the lower triangular matrices Q{z). Their partitioning brings us back to the two- 
diagonal pattern of eq. ( pT]) and replaces its scalars aj and Pj by the respective 
two-dimensional sub matrices Aj and Bj, 

I 



Q 



Ao 
So Ai 








5i ^2 



V 



In both the respective a = and a 7^ bases ([T0|) and (|T^ the Z^— plets of kets 
(|m -|- 1), |m -|- 2), . . . , |m -|- D)) with m = m{n) = nD — and with any may be 
denoted as In such an abbreviated notation our linear system @ implies the 

recurrence relations 



Fn = 



h. 



m{n)+l 



\ 



(^n) Bn-lFn-l, 



n 



1,2, 



(19) 



y hm(n)+D j 



which define the D— dimensional vectors of coefficients in terms of finite products 
of the certain D x D— dimensional matrices. In place of c/q = 1 in a consequently 
dimensional "democratic" partitioning we may use the shift do = D. Both these 



10 



options appear in our AS example where we recommended = D — p. The latter 
one is globally preferable as it leaves the uppermost element of Q vanishing, Aq = Q. 
The initial array Fq degenerates to the mere scalar norm then. 

At any the formal solution of the Schrodinger equation {H — E)\y) = 
may be re-written in the form of the double or partitioned sum, 

oo D oo 

b) = EE IH).[^n], = E (20) 

n=Oj=l n=0 

In a little bit vague sense it looks like an immediate hypergeometric-like general- 
ization of eq. (p!^. Equation (|19D defines all its coefficients in closed form. They 
depend on the "measure of asjTiimetry" a and on the unknown energy E = —k^. 



4 Example 

Our recipe strongly resembles the Hill- determinant method which proves useful in 



many (e.g., perturbative [^) applications. In the majority of similar applications 
one must analyze, first of all, the convergence of infinite series or (pO]). In 
X— representation their point-wise convergence is basically controlled by the asymp- 
totics of the coefficients. They are dominated by the purely kinetic terms which are 
asymptotically increasing. All the characteristics of the potential itself (e.g., parity 
mixing) will play, necessarily, a secondary role. 

The first non-trivial asymmetric potential \^['^'^l(x) seems best suited for a more 
explicit illustration of this role. Its coefficients hj = hj^\p) in both the p = and 
p = 1 solutions (^) are easily derived from the respective recurrences. Choosing the 
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simplest a = and using the same abbreviations aj and bj as above we have 



Q 



(0) 



U 








\ 


9 

u 
% 


ai 

9 «2 








n 




9 «4 











h 9 
-9 h 


as 

9 «6 















(21) 



Q 



(1) 










\ 


-5- ^1 


^2 









^3 9 

-9 h 












■■• / 



f22) 



V 

The contribution of the couphng g is clearly separated from the growing and energy- 
dependent kinetic terms. 

After a return to the general a 7^ we just have to modify the values of the 
matrix elements accordingly. We may preserve the reduction of bases (^) and ( ]TB| ) 
as well as their D = 2 partitioning. It is obvious that the exact Jacobi polynomial 
solutions may be reproduced in our D = 2 language. It is an instructive exercise to 
show how this reproduction proceeds. Firstly, the variability of the parameter a and 
of the energy or momentum k enables us to achieve a complete disappearance of the 
two-by-two submatrix Bk = at an arbitrary optional K. The resulting series (pO|) 
then strictly terminates and reproduces the known Gauss hypergeometric solution. 
Its termination just reflects the factorization of the secular determinant. 

Let us underline that the simpler, "termination-incompatible" basis (|10D with 
a = is an analogue of the non-WKB bases in Section 2.2. Hence, we may fix a = 
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and recall the same AS model also as one of the simplest illustrative examples of a 
general non-terminating solution. 

4.1 AS oscillator in the a = representation 

The AS solutions (||) may be split in the two separate sums with the well defined 
parity, 

|y[A5]^ _ |rl^^](j9)) = + \Y^"'^'^\p)). (23) 

The first few terms in the even partial sums with q = 1, 

|F(--)(0)) = ■ /i«(0) + IS5) ■ h^;\0) + ■ /.«(0) + . . . , 

)(!)) = IH3) • + ■ + \En) ■ hPil) + . . . , (24) 

as well as their odd, q = counterparts 

|r(°'^'^)(o)) = IH4) ■ hf\o) + iHg) • hf\o) + ■ hf\o) + ..., 

\Y^"''\l)) = • 4°^(1) + iSe) ■ /if (1) + |Sio) ■ /if (1) + . . . , (25) 
are easily computed in the recurrent manner, 

/if(0) = l, hf\0) = -g/a^, /if (0) = -feoM + //(aiaa), 

/if(l) = l, /if (1) = -^?M, /if (l) = -6i/a2 + ^?VK«2), .... (26) 

A compact general determinantal formula for these coefficients also exists [l^. It 
would enable us to re- write eq. (p3D, i.e., 

00 00 
\y^^'Kp)) = E |Sm.,p,i)) • h^^P) + E |S.a-i-i-P,P,o)) ■ /^f+i-p(p) (27) 

in the explicit form if needed. Here, we prefer the recurrent generation of the doublets 
of coefficients 



F'n+l-p — -^n+l-p(p) 



1^ ^2n+2-p(p) 



p = or 1, n = 0, 1, 
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as a matrix product, 

F,{p) = [~A,{p)]-' B.j{p) J = 1, 2, ... . (28) 

In our partitioned notation with D = 2 the solution \y) may be presented as a two- 
dimensional hypergeometric series since its matrix coefficients remain surprisingly 
elementary. 



1 I f 1 1 f Va2,+p-i 
l/a2j+p )[g 1 )[ 1 

As long as > ai > a2 > . . . at any k > 0, all our vectors of coefficients are well 
defined and unique. Their initialization is provided by the "model space" equation 
Ao{p)Fo{p) = which depends on p. Atp = 0, we have the vanishing scalar Ao{0) = 
while the exceptional singlet -fo(O) = h^o\o) (conveniently put equal to one) is the 
norm. In the parallel two-dimensional initialization aX p = 1, the ffist component 
4°^(1) = 1 of Fo(l) is the norm. The second component must be re-calculated, 
hl^\l)=gh^^\l)/i2K + l). 

We are ready to prove the convergence. Its decisive simplification occurs in the 
j ^ 1 asymptotic domain. The upper and lower components of eq. ( |28| ) decouple 
there in a p— independent manner. 



[F,^i{p)], , g = Oor 1. 

(29) 

For both our infinite series ( pT]) the proof is easy at x 7^ 0. As long as coshx > 1 for 
all the nonzero and real coordinates x, the ordinary geometric criterion together with 
the estimate (^) implies that our series ( P7D are convergent absolutely, i.e., for all 
the (complex) couplings g and energies — k^. The same geometric argument extends 
the validity of our conclusion to all the complex coordinates x + iy which lie out of 



a wiggly bounded domain such that | cosh(x + iy)\ = ysinh x + cos^y < 1 or, in a 
cruder approximation, out of the fairly narrow strip with |x| < ln(l + \/2) at least. 

On the real axis, an indeterminate behaviour of the type x cxd emerges at the 
point a; = 0. This follows from eq. (PUj) and from the slightly more sophisticated 
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Raabe criterion. Strictly speaking, this forces us to work on a punctured domain of 
X G (— oo, 0) U(0, c)o) in principle. As a consequence, logarithmic derivatives of our 
left and right Jost solutions have to be matched in the origin. This task is to be 
fulfilled numerically. Let us outline its two steps. 

4.2 Generalized parity 

Since our D = 2 hypergeometric AS series {x\Y^'^^^{p)) = ^^^\g, x, k) (pTl) satisfy the 
differential Schrodinger equation on a punctured domain (— oo, 0) 1J(0, oo) only, we 
necessarily have to match them in the origin. In the Poschl- Teller example of section 

2.3 where the non-matrix Gauss solutions also developed a certain discontinuity in 
the origin at a general unphysical energy E, the point has easily been settled after an 
account of parity. As long as our potentials lose their spatial symmetry in general, the 
parity is broken and a matching of the two sub-intervals (— oo, 0) 1J(0, oo) becomes 
nontrivial. 

We have to employ a broader invariance of our model(s) with respect to the 
product P of parity V with the reflections of couplings gj ~9j- The operator 
(such that = 1) commutes with our Hamiltonian(s), H = P H P. Each physical 
bound state ip{x) may be assigned an even or odd P— parity, Pipix) = ±'?/'(x). 

In a way resembling the parity-breaking systems with VT invariance the 
assignment of the P— parity to our AS states xio^ ^) depends on their normalization, 

{Px{9,x) = ±x{9,x)} =^ {P[9-x{9,x)] = t[9-x{9,x)]}. 

Fortunately, our AS coefficients h^n\p) = h[i\p^9) ^-^e explicitly defined by the 
triangularized Hamiltonians ( PT| ) and (^) and we immediately notice that 

hf{p. -g) = {-ir''^'hf{p,g). 

Both our AS hypergeometric-like series ip^^\g, x, k) = {x\Y^-^^\p)) (p^) behave as 
eigenstates of our double-parity operator P, 

P^^P\g,x,^) = ^(p){-g,-x,^) = {-ir^'~P\g,x,^) . 
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With a pair of some constants A4 M.{g) and M 7^ -^{g) we may postulate that 
the bound states read 

tlj^^^\x)= M^^^\g,x,K)+g-U^^^\g,x,K), x^Q. (30) 

The same (conventionally, even) P-parity may be assigned to all our physical solu- 
tions since their energy spectrum is non-degenerate. 

4.3 Match in the origin 

A return to the ordinary spatial parity V enables us to distinguish between the 
cosine-like (i.e., spatially even) and sine-like (i.e., spatially odd) components of our 
generalized hypergeometric functions (p^), 

c(a;, k) = i[(/?(°)(^,x, k) + (/?(°)(^,-x, k)], 

s{x,i^) = ^[f^^\g,x, k) - <^^^\g,-x, k)], 
c{x, t) = ^[v^^^Kg, X, k) + '^'■^\g, -X, k)], 
s{x,k) = ^[if^^\g,x, k) - <^^^\g,-x, k)]. 

The tildas ~ marking the asymptotical subdominance are not too relevant since we 
dwell in a vicinity of the origin where x = ±e ~ 0. Wavefunctions must be continuous 
there. 

The even, cosine-like components of our solutions satisfy such a requirement identi- 
cally. In the light of eq. (^) we are left with a reduced continuity condition 

Ms{e,K) + g ■Afs{e,K,) = 0, e ^ 0. (31) 

In the same manner, the continuity of derivatives is required. In the upper-case 
notation with abbreviations 

S{x,k) = ^[d^^^^\g,x,K) + d^<^^^\g,-x,K)], 
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C{x, 1^) = ^ [dx'p''^^ {g, X, k) - d^ifi'-^^ {g, -x, k)] , 
S {x, = ^ [dxv''^^ {g, X, k) + d^ip^^^ {g, -x, k)] , 
C{x,k) = ^[d^ip'^^\g,x,K) -d^ip^^\g, -x,k)] 
this leads to the second reduced matching condition 



MS{e,K) + g ■AfS{e,K) ^0, e ^ 0. 



(32) 



In the hmit £ — > a root K{e) of the two-dimensional secular equation 



det 



^ s{e,K,) s{s,K,) ^ 







\^ S(s,k) S{s,k) J 

will determine the physical energy. Matrix elements of this secular equation are 
convergent series in t — cosh~^ £ < 1, 

oo oo 



n=l 



n=0 



Sis, K)^J2i'^ + 2n) /i« (0, g) , §{8, = ^ (k + 2n + 1) /i« (1, ^) T 



n=0 



n=0 



Norms hQ\o,g) — hQ'\l,g) — 1 are fixed and the higher coefficients carry the 
dependence. An analogy with the spatially symmetric Poschl-Teller example of 
section 2.3 is fully restored. 
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5 Summary 



We described a new approach to the Schrodinger bound-state problem with any 
Rosen- Morse-hke multi-term potential (|l]). For all these forces we have shown how 

• the ordinary differential Schrodinger equation for the wave functions ^(x) may 
be reduced to a linear homogeneous algebraic problem Q{E)h = 0; 

• an "inspired" choice of the Lanczos-like (i.e., Hamiltonian-dependent) basis 
makes the related infinite-dimensional secular determinant vanish identically, 

det Q{E) = 0; 

• the very special (viz., lower-triangular) structure of our quasi-hamiltonian ma- 
trices Q{E) reduces the construction of the separate Taylor-like coefficients hn 
in our wave functions ipi^) to the mere (partitioned) two-term recurrences. 

On a characteristic AS example we have illustrated that 

• all our solutions ip{x) are convergent and may be understood as a certain 
generalization of the Gauss hypergeometric series (which further degenerates 
to the Jacobi polynomials at the physical energies in the solvable cases); 

• a certain generalized parity symmetry of our forces enables us to determine Jost 
solutions which are compatible with both our asymptotic boundary conditions; 

• via our final two- by-two condition (|31|) + (|32D, the values of the remaining two 
free parameters (viz., energy and p— mixing) in our Jost solutions may (and 
have to) be tuned to their necessary continuity and smoothness in the origin. 
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Table 1: Shape-invariant potentials on the line 



model 


V{x) 


V{-oo) 


V{oo) 


polynomial ipi^) 


harmonic 




oo 


oo 


Laguerre 


Morse 


ae~^ + 6e~^^ 


oo 





Laguerre 


Rosen-Morse II 


/ / cosh^ X + g tanh x 


-g 


g 


Jacobi 


scarf II 


{f + g sinh x) / cosh^ x 








Jacobi 
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